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Context. The combination of image restoration and a Fabry-Perot interferometer (FPI) based instrument in solar observations results 

in specific calibration issues. FPIs generally show variations over the field-of-view, while in the image restoration process, the 1-to-l 

relation between pixel space and image space is lost, thus complicating any correcting for such variations. 

Aims. We develop a data reduction method that takes these issues into account and minimizes the resulting errors. 

Methods. By accounting for the time variations in the telescope's Mueller matrix and using separate calibration data optimized for 

the wavefront sensing in the MOMFBD image restoration process and for the final deconvolution of the data, we have removed most 

of the calibration artifacts from the resulting data. 

Results. Using this method to reduce full Stokes data from CRISP at the SST, we find that it drastically reduces the instrumental 
and image restoration artifacts resulting from cavity errors, reflectivity variations, and the polarization dependence of flatfields. The 
results allow for useful scientific interpretation. Inversions of restored data from the 6 sunspot AR11029 using the Nicole inversion 
code, reveal strong (~10 km s~^) downflows near the disk center side of the umbra. 

Conclusions. The use of image restoration in combination with an FPI-based instrument leads to complications in the calibrations 
and intrinsic limitations to the accuracy that can be achieved. We find that for CRISP, the resulting errors can be kept mostly below the 
polarimetric accuracy of ~10~^. Similar instruments aiming for higher polarimetric and high spectroscopic accuracy, will, however, 
need to take these problems into account. 
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1. Introduction 

As technical possibilities increase, more advanced and compli- 
cated instruments are being designed for observational studies 
of the Sun. New instrumentation is generally designed to offer 
better spatial resolution, better spectral coverage, and full polar- 
ization information. Higher spatial resolution will help resolve, 
e.g., small-scale magnetic fields in the quiet Sun. In sunspots, the 
magnetic fields and flows in penumbra, lightbridges, and umbral 
dots make interesting targets. Better spectral coverage and polar- 
ization information are needed to enable meaningful interpreta- 
tion of the observed features. Adaptive optics (AO) are currently 
available at the major ground-based solar observatories and have 
greatly improved the performance of telescopes. 

Additionally, postprocessing of high-cadence imaging data 
to correct the blurring eflTect of seeing has made a huge impact 
on image quality. Image restoration methods such as MOMFBD 
(Multi-Object Multi-Frame Blind Deconvolution, Ivan Noort 
et al.|2QQ5}|Lofd ahl 2002) and Speckle interferometry (see 



In this context, Fabry-Perot interferometer (FPI) based in- 
struments have become quite popular a nd are in operation at 
many solar observatories, e.g. CRISP (Scharmer et al.|[2008|) 



von 



der Luhe||1993| and references therein) have become essential 



for ground-based, high spatial-resolution solar imaging. Because 
these methods allow combining all the images obtained within a 
dataset, while retaining high image quality, it removes the need 
for frame selection and results in more reliable and higher signal- 
to-noise data. 
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at the Swedish 1-m Solar Telescope (SST), the Telecentric 
etalon solar spectrometer (TESOS, |Kentischer et al.|1998), and 
Gottingen Fabry-Perot spectropolarimeter ( Puschmann et al. 
2006 ) at the Vacuum Tower Telescope and the Interferometric 
bidirectional spectrometer (IBIS, [ Cavallini||2QQ6 ) at the Dunn 
Solar Telescope. Similar instruments are likely to be developed 
for the future 4-m class solar telescopes ATST and EST. The 
main advantages of FPIs over Lyot filters are that they have a 
high transmission, allow for rapid wavelength tuning, and allow 
for a dual-beam polarimetric setup with a polarizing beam split- 
ter close to the final focal plane. A disadvantage is that even 
very well-made FPIs will show variations in the width and cen- 
tral wavelength of the transmission peak over the field of view 
(FOV, e.g. [Reardon & Cav allini'^8). When the FPIs are set 
up in a telecentric mount, as with CRISP, these variations cause 
intensity variations on relatively small spatial scales. In colli- 
mated mount, such small-scale features disappear, but the central 
wavelength shifts systematically away from the center of the im- 
age. For further comments on FPIs in telecentric vs. collimated 
mount see Scharmer^ (2006 J . 

We plan to discuss a theoretical approach to the processing 
of full- Stokes line scans from FPI-based instruments in a future 
paper, but here we present the practical application of that ap- 
proach to real CRISP data of a sunspot observed in the 6302.5 A 
Fe I line, including all the necessary data reduction steps, such as 
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correction for telescope polarization, flatfielding, image restora- 
tion, and demodulation. In Sect. [2] we discuss the general setup 
of the data acquisition procedures, and in Sect. [3] the challenges 
and possibilities for calibration procedures. The proposed cali- 
bration procedure is presented in Sect. [4] results of the procedure 
applied to CRISP data are presented in Sect. [5] and the conclu- 
sions in Sect. [6l 



2. Polarimetric data acquisition 

2.1. The optical setup 

The data acquisition setup used wit h CRISP at the SST (see 
Fig. \]) is similar to the one used by |van Noort & Rouppe van 
derVoort| ( [20Q8] ), with the notable difference that the SOUP Lyot 
tunable filte r has bee n replaced with the CRISP Fabry-Perot fil- 
ter fScharmer et al. 2008 ). The significantly higher (~ 8x) peak 
transmission of the Fabry-Perot filter results in a higher S/N for 
similar exposure times than with the Lyot filter, and the use of 
a dual beam polarimetric setup results in a significant reduction 
of the seeing-induced crosstalk. In addition, the higher tuning 
speed of CRISP allows for observing programs scanning over 
-10 spectral points instead of the two to four points one would 
typically get with the Lyot filter. A downside, however, is that 
the placement of the filter cavities in a telecentric mount con- 
figuration leads to small-scale pass-band shifts and transmission 
profile width variations across the FOV, due to cavity and reflec- 
tivity errors. In the case of CRISP, the high- and low-resolution 
FPIs have been placed ~10 cm away from the focal plane in or- 
der to smooth out such variations. 

CRISP data are obtained while continuously cycling through 
a 4-state liquid-crystal (LC) scheme, which converts combina- 
tions of the incoming polarization states Q, U, and V to linear po- 
larization that can be analyzed with the polarizing beam splitter 
close to the final focal plane. Data are simultaneously recorded 
using two narrow-band cameras (called transmitted and reflected 
camera) in a dual-beam setup and a wide-band camera (receiv- 
ing 10% of the light passing the prefilter via a beam-splitter) 
placed before the liquid-crystal variable retarders. The frame rate 
of -36 Hz is set by the rotating chopper, which has a duty cycle 
setting the exposure time at 16 ms. CCD readout is performed 
during the dark part of the cycle. For the data shown in this pa- 
per, 4 full LC-cycles were completed before tuning to the next of 
a total of 12 wavelengths. The complete linescan of 4 LC-states 
X 4 repeats x 12 wavelengths = 202 frames covers ~ 6 seconds. 



2.2. Calibration data 

The data required for calibrating the science data from CRISP 
are obtained in the following way. The dark images are recorded 
in a similar way as the science data, but with the light path 
blocked after the exit window of the vacuum tube. The final dark 
is the average of a large number of darks (typically 100-1000). 
Flatfields are obtained by summing many (-1000) images of the 
Sun taken with the telescope pointing moving across the solar 
disk. As such, the flatfields include both transmission variations 
over the FOV as actual CCD gain variations, as do the gains that 
are derived from these flatfields. Therefore, in this paper, "flat- 
fielding" and "gain correction" imply the same correction. The 
polarization calibration data is obtained by producing light with 
known polarization states using a linear polarizer and a quarter- 
wave plate placed below the exit window of the vacuum tube in 
the tower. The modulation scheme is then calibrated using the 
method described in Appendix |A| The (time-dependent) tele- 
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Fig. 1. Schematic drawing of the optical set-up at the SST. After 
reflecting ofl' the dichroic mirror separating the red and blue 
beams (not shown), the red beam passes the chopper and pre- 
filter. Then a beam splitter sends 10% of the light to the WB cam- 
era. The other beam passes the LCs, the high- and low-resolution 
etalons of the CRISP filter (HRE and LRE resp.), and finally a 
polarizing beam splitter separating the p- and s-polarized light 
to the transmitted and reflected beams. All three cameras are 
1024x1024 pixel CCD (Sarnofl) cameras running at -36 Hz. 



scope Mueller matrix has been measured by Selbing ( 2006| ) and 
is assumed to be known here. 



2.3. Image restoration 

For ground-based, high-resolution imaging, image restoration 
is a vital tool. However, it is especially important for multi- 
wavelength and polarimetric observations, because it will reduce 
difl'erential image shifts and blurring between images recorded 
at difl'erent wavelengths and polarization states. For polarimetric 
data, this implies that the image restoration has to be done be- 
fore demodulating the data (to convert the LC-states to Stokes 
parameters I,Q,U,V), because the images of the difl'erent LC- 
states are each associated with a particular seeing disturbance 
that must be compensated for individually, to reduce artificial, 
seeing-induced signals. In o ur analysis we focused on image 
restoration using MOMFBD (van Noort et al. 2005 ), but some 
of the complications discussed below are equally relevant for 
Speckle restorations. 

MOMFBD has several properties that can complicate cali- 
brations. 



1. The direct relation between an image pixel in the original 
data and the restored image is lost. 

2. Artificial features in images can influence the wavefront 
sensing and can be amplified by the deconvolution. 

3. Input images are assumed to be monochromatic, although 
MOMFBD can handle diff'erent wavelength channels within 
one restoration. 

4. Noise is assumed to follow a Gaussian distribution with a 
constant width over each isoplanatic patch, which can be a 
problem when high-contrast targets are observed. 

The loss of the direct relation between input data pixels and 
restored data pixels means that all calibrations that cannot be 
done before image restoration, and hence before demodulating, 
cannot be done on a pixel-by-pixel basis, but involve the point- 
spread-functions (PSFs, which result from seeing and the finite 



Schnerr et al.: Stokes imaging polarimetry using image restoration: A calibration strategy for Fabry-Perot based instruments 



3 



pupil size) of the input images. Similarly, artificial features in 
the input images should be removed as much as possible before 
image restoration, because they will interfere with the wave- 
front sensing and can be strongly amplified by the deconvolu- 
tion. Input images that are not monochromatic over the FOV are 
rather hard to deal with, because this cannot be properly cor- 
rected without knowing the entire line profile. 

Both the wavelength changes over the FOV and the assump- 
tion of noise with a Gaussian distribution of a constant width are 
limiting factors for the maximum (polarimetric) accuracy that 
can be obtained when using image restoration. To avoid interfer- 
ence from artifacts in the narrow band cameras on the wavefront 
sensing, one can only use images from the wide band camera 
for the wavefront sensing and apply the corrections to the data 
of all cameras. However, the additional wavefront information 
contained in the narrow band data is then not used for the image 
restoration. 

The changing beam configuration due to seeing corrections 
by the AO (and tip-tilt mirror) in principle also aff'ects data cali- 
bration, because it changes the way the optics in front of the AO 
are sampled. However, it does not prevent pixel-based calibra- 
tions, and since we assume the telescope matrix is constant over 
the FOV, the eff'ects of the changing beam configuration are not 
considered here. Image restoration of the data used in this paper 
was performed with MOMFBD assuming isoplanatic patches of 
96x96 pixels (-5.7x5. 7'') and using 37 Karhunen-Loeve modes 
for modeling the wavefront phases. 

3. Instrument calibration 

A straightforward calibration scheme typically involves the fol- 
lowing steps. 

1. Dark-correct data from the diff'erent cameras (transmitted, 
reflected, and wide band). 

2. Gain-correct the difl'erent LC- state and wavelength images 
of all cameras. 

3. Use image restoration to combine the data for the difl'erent 
cameras and states, and correct for the seeing. 

4. Demodulate transmitted and reflected beams using the polar- 
ization calibration matrix (see Appendix [A). 

5. Combine Stokes images from the transmitted and reflected 
beams to reduce seeing crosstalk. 

6. Compensate for telescope polarization using the inverse of 
the telescope Mueller matrix. 

Due to specific instrumental properties, such as polarization- 
dependent flatfields and FPI inhomogeneities, polarimetric and 
spectral calibration problems arise in steps 2 to 4. These are dis- 
cussed in more detail below. 

3.1. Polarimetric issues 

For polarimetric calibration of the data it is convenient to split 
the calibration in two parts. 

1. Telescope, from the 1-m primary lens at the entrance of the 
vacuum tube down to the exit window at the bottom of the 
vacuum tube. Due to the 45° angle of incidence on the rotat- 
ing azimuth- and elevation-mirrors, this part is intrinsically 
time variable, but is assumed to be independent of the line- 
of- sight (constant over the FOV). A model describing the 
Mueller matrix of the telescope as a function of the azimuth 
and el evation angles at 6302 A was determined by |Selbing| 
( |2QQ6l ). 



2. Instrument, from the vacuum exit window to the CCDs. 
This part is FOV-dependent (e.g. dust close to the focal 
plane), but can be calibrated daily using calibration optics 
and is assumed not to vary in time. It is described by a mod- 
ulation matrix, which converts the input Stokes vector into 
intensities measured with the LC- states. How this matrix can 
be determined is described in Appendix [A[ 

We can describe the total conversion matrix of the system (Ji) 
as a matrix product of the Mueller matrix of the telescope T and 
the modulation matrix of the instrument Ai: 

J[(x,yj) = M(x,y)'T(t), (1) 

where the transmitted and the reflected beams have their own 
(although related) modulation matrices. Because the intensity of 
the light during the determination of the telescope matrix and the 
modulation matrix is not known, an additional gain calibration is 
still required. The telescope matrix T(t) varies with time, so the 
total system response matrix ^ is time dependent. This means 
that unless the flatfield images and the data are recorded at the 
same time, which is not practical, an appropriate gain correction 
(for the modulated data) must be computed. To this end, the flat- 
field must first be transformed to a time-independent location. 
The obvious (and only) place available for this is the plane of 
the sky, i.e. on the Sun, since we assume that all time variability 
is contained in the telescope matrix. The details of this procedure 
are discussed in Sect. lU 

3.2. Spectral issues 

An important calibration problem is that part of the intensity 
changes over the field of view are not the result of real CCD gain 
variations, but of variations in FPI reflectivity and cavity size 
(see, e.g., |Reardon & Cavallini|20Q8|[Cauzzi et al.|20Q9 ). When 
observing at a wavelength where SI/SA 0, e.g. in the wing of 
a spectral line, a wavelength shift due to a cavity error results in 
a change in the observed intensity that is completely unrelated 
to the gain of that particular pixel. Similarly, variations in the 
reflectivity cause changes in the integrated transmission profile. 

The question, however, is which flatfield correction should 
be applied to the data. Although the intensity variations due to 
cavity errors are not real gain changes, the same cavity errors 
also aff'ect the science data. The line profiles in the science data 
can be very diff'erent, though, so the intensity changes will in 
general not be the same as in the flatfield data. A change in wave- 
length scale due to a cavity error just cannot be corrected in a 
flatfielding procedure, which only changes the intensity, unless 
one knows the exact line profile. For a data calibration scheme 
without image restoration, one could attempt to somehow de- 
termine the real pixel gains, apply those, and account for the 
wavelength shifts later on in the analysis. For image restoration, 
however, the intensity variations that then remain could pose a 
problem for the wavefront sensing and/or the deconvolution. 

The strategy that we propose is discussed in the next section. 

4. Calibration scheme 

For the wavefront sensing part of the image restoration, it is im- 
portant to use images that are as monochromatic as possible and 
that have no sharp artifacts. Sharp features can result from CCD 
pixel rows that have diff'erent responses to the polarization of the 
incoming Stokes vector, when the S/P beamsplitter (see Fig. [T]) 
is not 100% eff'ective. These pixels then stand out in some of 
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the matrix elements of M. Such features should be reproduced 
in flatfields that are properly converted to the time of the obser- 
vations, because this information is contained in the modulation 
matrix, as long as the polarization of the incoming light is low, 
the telescope polarization is known, and the instrument is sta- 
ble. Therefore these features should be removed from the data 
by flatfielding. 

Dealing with the non monochromatic nature of the data, re- 
sulting from FPI transmission variations over the FOV, is more 
complicated. In the final restored data, one can attempt to take 
the wavelength shifts caused by cavity errors into account, al- 
though as a result of the image restoration process the eff'ective 
shift in a given pixel in the restored image is always some mix- 
ture of the shifts of all CCD pixels within the PSFs of the dataset. 
For the image restoration, the intensity changes are more prob- 
lematic, especially if they contain high spatial frequencies, be- 
cause these distort the wavefront sensing. Intensity changes due 
to reflectivity and cavity errors of the FPI cannot be interpreted 
as simply brighter or darker regions in the observed object by 
the wavefront sensing, because they do not move with the object 
while the seeing varies throughout a dataset. 

Our approach is to use flatfields that make the data as smooth 
as possible, before using the data for the wa vefro nt sensing. This 
is done using the recipe described in Sect. |4.1[ These flatfields 
include intensity changes due to varying cavity properties over 
the FOV. This is a straightforward choice, but a relatively ar- 
bitrary one, because such a flatfield will only properly correct 
these intensity changes if the observed line profile is the same as 
the one of the flatfield data, which is in general not the case. For 
the final restored data however, we do not want to use these flat- 
fields, because they modify the observed line profiles, but rather 
those flatfields that are properly corrected for the varying cavity 
properties. 

A relatively simple way to obtain flats that are not afl'ected 
by cavity errors in the FPI would be to use a wavelength close 
to the observed line where 61 ISA = 0. However, such a method 
does not account for the wavelength dependence of the flatfields 
from variations in the prefilter transmission profile, LRE/HRE 
co-tuning errors and the wavelength dependent fringes that are 
observed in the Sarnoff' cameras. Additionally, the wavelength 
where such flats have to be taken might not correspond to the 
preferred continuum point, because the instrumental response 
has to be taken into account, increasing the time required for 
calibrations. Therefore we developed a more accurate routine to 
derive proper flatfields, which simultaneously solves for the FPI 
cavity and reflectivity errors, and prefilter changes over the FOV. 
This routine is described in Sect. 14.21 



4.1. Time-independent and image-restoration flatfields 

The proposed scheme thus consists of creating flats especially 
designed to minimize image restoration artifacts and time- 
independent flatfields that should properly flatfield the data with- 
out modifying the line profile. The image restoration flatfields 
are applied before the image restoration process, but removed 
again afterwards when the time-independent flatfields are ap- 
plied. 

To obtain time-independent flatfields, we have to invert the 
modulation matrix of each pixel to obtain the demodulation ma- 
trix (Al"\ the inverse of the modulation matrix At). By apply- 
ing the demodulation matrix to the flatfield data of the 4 LC- 
states, we obtain flatfield images for all Stokes parameters under 
the vacuum-tube's exit window. Using the inverse Mueller ma- 
trix of the telescope, valid for the time at which the flatfields 



were taken, we can calculate the Stokes flats on the Sun (on the 
plane of the sky). The flatfields represent the eff'ective flatfields 
for Stokes I, Q, U, and V. An overview of how both the time- 
independent and image-restoration flatfields are determined is 
shown in Fig. [2] 

To calculate the flatfields relevant for the time the observa- 
tions were obtained, one has to apply the telescope matrix ap- 
plicable for that time to the flatfields on the Sun and then ap- 
ply the modulation matrix. For this we only use the Stokes I 
flatfield (and set Stokes Q, U, and V to 0) to avoid imprinting 
a polarization pattern on our LCI.. 4 flats (which results from 
the FOV- variations in the Stokes QUV-to-LC1..4 components of 
At). These flats will include intensity changes due to cavity and 
reflectivity errors, which should smooth these variations in the 
data (when wavelength shifts due to line-of- sight velocity vari- 
ations are small compared to the cavity errors). Also the polar- 
ization levels should be properly reproduced to the accuracy that 
the telescope Mueller matrix and modulation matrix are known, 
if the intrinsic polarization of the observed target is low. 

4.2. Determining ttie proper flatfields 

The FPI has small variations in the cavity size and reflectivity, 
which results in changes in the central wavelength and width 
of the transmission profile, respectively. In addition, the central 
wavelength of the prefilter shows a slight variation over the FOV 
as the angle of incidence on the prefilter changes (due to the 
physical dimensions of the CCD). To determine the real trans- 
mission over the FOV, one has to account for these eff'ects. At 
the SST, flatfields are determined by summing many (-1000) 
exposures per state and wavelength obtained while the tele- 
scope pointing is moving across the solar disk (to avoid intensity 
changes over the FOV due to limb darkening). If one assumes 
that the observed line profile of the averaged flatfield is constant 
over the FOV, then the diff'erent factors can be determined by 
using an iterative fitting scheme. 

The transmission of a single FPI is given by the Airy function 
(e.g. |Kentischer et al.|1998l ): 



1 



f FPI 



1 -\- Fsm^(2ndncose/A) 



where the finesse F is given by 



4R 



(l-7?)2 



(2) 



(3) 



with R the reflectivity, d the cavity size, n the refractive index 
in the cavity, the angle of incidence, and A the wavelength. 
The total transmission ^ of the dual FPI setup of CRISP is cal- 
culated by multiplying the transmission of the high-resolution 
(Jhr) and low-resolution FPI (Ti^), assuming ^ = and assum- 
ing that cavity and reflectivity errors in the low-resolution FPI 
can be neglected because those in the high-resolution FPI will 
be dominant: 



^^(R,AR) = T,,(R)'TUR^m- 



(4) 



Because the spectral range used is only a small fraction of the 
prefilter transmission profile width, we can limit the number of 
free parameters by not including the average prefilter shape in 
our fits, but instead only fit the deviations from the mean prefilter 
shape, assuming a linear dependence on wavelength. 

Given the constant line profile (p(A), average transmission 
profile ^(R) of the combined FPIs, and (wavelength indepen- 
dent) gain g, the observed line profile in pixel / will be given 
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Fig. 2. Flowchart illustrating the procedure used to calculate the flatfields used for wavefront sensing and the time-independent 
flatfields corrected for wavelength shifts. 



by 

^ohsA^) = [(1 + ^p/,/) • W - ^o)] • ^(/? + A/?,) * - A^,) (5) 

where gi is the pixel gain, D^f^i the prefilter deviation factor, 
AT?/ the reflectivity error, A/1/ the cavity error, and the central 
wavelength of the line. 

The iterative scheme is initiated by setting gi to the inverse 
of the mean intensity of the line profile, - A/1/ to the mini- 
mum of a parabolic fit to the line minimum, and Dpf^i and A/?/ to 
zero. After applying these corrections, one can overplot all the 
observed line profiles and determine the average line profile by 
calculating a spline through a set of points at fixed wavelengths, 
for which the intensity is determined by a;^^ minimalization rou- 
tine (see Fig. [3} all fitting was done using the mpfit package, see 
Markwardt|2QQ9l |More| 1978] ) . The "real" line profile is then cal- 
culated by deconvolving this fitted line profile with the average 
transmission profile. As we now have an estimate for the aver- 
age line profile and transmission profile, we can fit for the four 
unknown parameters in Eq. |5]for each pixel. These parameters 
allow us to improve our estimate for the average line profile and 
transmission profile. These two steps can be repeated until no 
more improvement is made. 

Results of this fitting routine applied to a real dataset (see 
Sect. [5]) are shown in Figs. [3] and |4] The total gain of a pixel 
is determined by the wavelength-independent part gi and a 
wavelength-dependent part given by the ratio of the observed 
intensity and the fitted line profile, as shown in Fig. [3] 



4.3. The prefilter transmission profile 

After the prefilter deviation factor and average line profile have 
been determined, we use the average line profile of all the pix- 
els (corrected for cavity errors and gain) to determine the av- 
erage prefilter shape Tpp. The prefilter, which is mounted at a 
slight angle to tune the central wavelength, is assumed to have a 
Lorentzian profile: 



7^pf(^) = 



1 



1 + (2(/1-/1o,pf)/w> 



2A^ca 



(6) 



where /1o,pf is the prefilter central wavelength, w the full- width 
at half-maximum (FWHM), and A^cav(=2) the number of cavities 
of the prefilter. 




Fig. 3. Density plot of intensity (arbitrary units) vs. wavelength 
relative to line core for all points in the summed flatfield data 
before the first iteration (top, where the line center has been de- 
termined using a parabolic fit to the line core) and the same plot 
after the second iteration (bottom). The average line profile (f)(A), 
described by a spline curve through several nodes (red plusses), 
is overplotted in red. 



The prefilter parameters /1o,pf and w are determined by 
matching the average observed line profile to the line profile 
of the FTS atlas (the atlas acquired with the Fourier transform 
spectrometer at the McMath-Pierce Telescope, available from 
Brault & Neckel|1987 ) convolved with the average transmission 
of CRISP (see Fig. |5]), masking the points in the telluric line. We 
find a prefilter width of 4.7+0.4 A, which agrees with the speci- 
fication of 4.4-4.6 A, and a central wavelength of 6301.2+0.2 A. 
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Fig. 4. Wavelength shift resulting from cavity errors (left), variations in the FPI reflec tivit y (middle, nominal value 93.5%), and the 
prefilter variation factor (right) as determined by the fitting routine described in Sect. 42^ (see also Eq.[5]) when applied to the 6302 
A dataset presented in Sect, [s] 




Fig. 5. The prefilter shape (dashed line) is determined by match- 
ing the uncorrected average spectrum (crosses) to the FTS con- 
volved with the CRISP transmission (full black line). After cor- 
recting for the prefilter, an excellent match is found (red circles). 
The points in the telluric line (gray circles) were not used in the 
procedure. 

4.4. Iterative versus direct method 

After calculating both the actual CCD gains that we want to ap- 
ply in the end and flatfields that should be used for the image 
restoration, we still have to account for the impact of the image 
restoration process, as the result of which the 1-to-l pixel re- 
lation is lost. Here we discuss two possible approaches to that 
problem: 

- Direct method. Apply the flatfields prepared for the image 
restoration to the data and after the image restoration remove 
these flatfields and apply the time-independent flatfields, tak- 
ing the average of the PSFs determined by the image restora- 
tion process into account. 

- Iterative method. Use the PSFs that have been determined 
by the image restoration to iteratively determine the Stokes 
images that best reproduce the observed data given the tele- 
scope Mueller matrix, modulation matrix, PSFs, and gains. 

An overview of both methods is presented in Fig. [6] For the 
direct method, flatfields for the time of the observations are de- 



termined as described in Sect. 4.1 which are used for the image 
restoration process. These flatfields are then removed from the 
restored data and the true flatfield applied by multiplying the 
data with the ratio of the image restoration flatfield and the true 
flatfield, convolved with the average PSF. The demodulation is 
then performed with the inverse of the modulation matrix con- 
volved with the average PSF. The idea behind this is that each 
pixel in the restored image "samples" the whole PSF, and there- 
fore the convolved calibration data should be applied. 

An alternative method, which is simpler but computation- 
ally more expensive, is the iterative method, which will be de- 
scribed in more detail in a future paper we are preparing. This 
approach solves the entire forward problem, including all the 
calibrations and individual PSFs, and iteratively improves the in- 
coming Stokes images to find the solution that fits the observed 
data best. 



5. Results 

5.1. Example: sunspot 6302 

As a test of our reduction method, we reduced a dataset from 27 
October 2009, which has part of the S sunspot AR11029 in the 
field of view (see Fig. [7]) at a heliocentric angle of about 45° (ji ^ 
0.7). The data contain 11 equidistant wavelength points around 
the core of the Fei 6302.5 A line with a wavelength separation 
of 48 mA and a continuum point at -1-600 mA. 



5.2. Stokes parameters 

In Fig.|7]we show the improvement in the Stokes parameters that 
is reached by using the scheme described in Sect.|4]as compared 
to straightforwardly applying standard calibrations as listed in 
Sect. |3] Polarimetric errors due to the changing telescope polar- 
ization that influences the flatfields are visible as saddle-shaped 
off'sets, a vertical fringe-like pattern, and sharp lines in Stokes 
Q, U, and V. These errors almost completely disappear with the 
new calibration method. 
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Fig. 6. Flowchart illustrating the direct and indirect method to calculate the final Stokes I, Q, U, and V images using image restora- 
tion. 



The origin of these problems lies in the variations in the mod- 
ulation matrix. In Fig. [8] we show two representative examples of 
the 16 components of the modulation matrix: the contributions 
of Stokes I and Q to LC-state 1. As for all the LC-states, the 
Stokes I contribution shows a vertical fringe-like pattern over the 
whole FOV, some smaller scale patterns, likely related to dust, 
and sharp lines. The Stokes Q contribution is much smoother 
and mainly shows a saddle shape and some sharp lines, which 
is typical of the contributions of Stokes Q, U, and V to the LC- 
states. Because such patterns are present in the modulation and 
demodulation matrices, it is not possible to remove them both 
from the Stokes I, Q, U, and V images and from the correspond- 
ing modulated LC-states at the same time. Flats being chosen to 
make the Stokes images smooth, implies that the modulated im- 
ages will not be smooth, and vice versa. Therefore our scheme 
uses two sets of flats, one to make the modulated images smooth 
for the wavefront sensing and one to properly flatfield the Stokes 
images. 

The dark regions in the Stokes V map in Fig. [7] are actual 
signals redshifted into the continuum wavelength due to strong 
downflows (see also Sect. 5.3). The small spot of fringes near 
the bottom of the Stokes Q, U, and V images is likely from a po- 
larizing speck, e.g. dust, in the telescope, which is not accounted 
for in the telescope model. Outside these areas, the noise in the 
continuum point of the Stokes I normalized polarization images 
is -0.3%. 



5.3. Inversion results 



The data are inverte d using the inversion code Nicole ( |Socas- 
[Nav arro et al.| [20TT] ). For the inversions we use three nodes in 
temperature, one node (i.e. the parameter has a constant value) 
in line-of- sight velocity (vios), and one node for the line-of-sight 
magnetic field strength and the microturbulence. The average 



transmission profile (see Sect. 4.2) is used to account for the 
finite spectral resolution of CRISP. 

Results of the inversions are shown in Fig. [9] The shift in 
the velocities due to cavity errors of the FPI is corrected by sub- 
tracting a map of the cavity errors as determined in Sect. 4.2 
convolved with the average PSF of the whole dataset. The zero 
velocity reference was set to the average velocity in the umbra. 
Adopting this velocity reference, we find a convective blueshift 
as measured from the granulation in the FOV of 0.3 kms"\ 
which is consistent with the theoretical value of 0.4 kms"^ (de 
La Cruz Rodriguez et al.|20lT ). 

Both the velocity and longitudinal magnetic field can reliably 
be determined from the data, showing fast (~5 kms"^) outflows 
at the center-side penumbra and downflows (up to -10 kms"^) 
on the intersection between the umbra and penumbra slightly 
below the middle of the FOV. These downflows will be discussed 
in more detail in a future paper. 



6. Conclusions 

We have presented a method for reducing full Stokes data ob- 
tained with CRISP at the SST, an instrument with a fast- tuning 
Fabry-Perot filter in a telecentric mount. Although some calibra- 
tion issues are inherent to such an instrument used together with 
image restoration, most problems can be overcome by using the 
proper calibration procedures. 

The impact on the flatfields of wavelength shifts over the 
FOV, owing to cavity errors, can be corrected for. This is impor- 
tant for preventing changes in the line shape caused by errors 
in the gain tables. For instruments with the FPI in a collimated 
mount, the cavity errors occur in pupil space, and are therefore 
not a problem in this respect. These instruments do show radial 
wavelength shifts over t he FOV that have to be accounted for 
( [Janssen & Cauzzi|2006| ). 
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Fig. 7. Continuum Stokes I, Q, U, and V images (from top to bottom) of the sunspot observed with the SST in the 6302.5 A Fei 
line on 27 October 2009. The images on the left were calibrated with the described calibration scheme and the iterative method, 
the images in the middle used the direct method instead, and the images on the right were calibrated by straightforwardly applying 
standard calibrations as listed in Sect.[3]and using the direct method. The full range of the polarized images is -1.5% to +1.5%. The 
arrows indicate where sharp artifacts occur for some methods. Continuum offsets in the Q, U, and V images, most likely resulting 
from discrepancies in the telescope model, of up to -0.35% have been subtracted. The small dark region visible in the Stokes V 
image is actual signal in a region with strong downflows (see text). 
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Fig. 8. Variations in modulation matrix elements over the FOV. The element describing the contribution of Stokes I to LC- state 1 is 
shown on the left and the contribution of Stokes Q to LC- state 1 is shown on the right. 
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Fig. 9. Results of inversions of the data, showing normalized continuum intensity {left), the line-of-sight component of the magnetic 
field (middle), and the line-of-sight component of the velocity {right). The inversions reveal quite strong outflows (up to about 5 
km s"^) and strong downflows (~10 km s"^). The limb is towards the upper right. 



Without a proper calibration scheme, the polarimetric im- 
print of the flatfielding can cause strong saddle- shaped off'sets, 
fringes, and sharp lines from deviating pixel rows/columns that 
are enhanced by the deconvolution in the image restoration. By 
using the proposed calibration scheme these artifacts can almost 
be completely removed. Any remaining weak artifacts are most 



likely due to discrepancies in the SST telescope model, which 
is only accurate to a few tenths of a percent, and nonzero in- 
trinsic polarization. The severity of the artifacts depends on the 
(variation in the) telescope polarization and the variations in the 
modulation matrix elements over the FOV. For telescopes that 
have low and/or constant intrinsic polarization and instrumenta- 
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tion whose modulation matrix elements are nearly constant over 
the FOV, such flatfielding problems will be much less severe. 

Velocities, as determined from inversions of the data, can be 
corrected for cavity errors using a convolved version of the cav- 
ity map. However, the PSFs will vary from wavelength to wave- 
length resulting in a non equidistant wavelength sampling. Both 
the wavelength sampling and the variations in the width of the 
FPI transmission peak will vary from pixel to pixel, but can in 
principle be taken into account during the inversions. An intrin- 
sic limitation is the use of non monochromatic input images for 
the image restoration, which results in deconvolution artifacts: 
mixing of information from different wavelengths through the 
image restoration procedure. 

Modem developments in FPI instrumentation, such as rapid 
(kHz) modulation combined with charge shuffling (which may 
allow for demodulation before image restoration), could reduce 
calibration problems. However, to reach the polarimetric and 
spectroscopic accuracy that the future large solar observatories 
EST and ATST are aiming for with the next generation of Fabry- 
Perot systems, advanced calibration techniques and optimized 
optical designs to deal with the issues discussed here will be re- 
quired. 
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Appendix A: Determining the modulation matrix 

When fitting the pixel-dependent system response matrix At, it is important to reaUze that the incoming intensity / is not known. A fit can therefore only be made 
on the fractional polarization, which means that we need to normalize the data with the intensity. The resulting data, lacking absolute transmission information, can 
only be used to determine the modulation matrix up to a constant factor A~^. Moreover, the normalization of the data requires us to compute the intensity from the 
data using the inverse of the unknown modulation matrix At. It is therefore more natural to formulate the fit directly in terms of the demodulation matrix Ai~^ : 



M- 
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The calibration data is taken using a linear polarizer followed by a /1/4-wave plate, which produces a known (relative) polarization, independent of polarization 
of the incoming Hght, by separately rotating the two plates. Given the calibration data v and the known (normalized) Stokes vector of the incoming light S, we have 
for each calibration measurement point n with unknown intensity 



bijVj,n = S inh = S i^n ^ bikV^^n 



(Al) 



so that for a given guess-solution for At ^ the total accumulated fit error between the data and the guess matrix can be written as 



A ( A A \ 

i=\ n \j=\ j=\ ) 



(A.2) 



Due to the use of normalized data, we have 



Substitution in Eq. 



A.2 



E 



bii = A. 



eliminates the contribution to;^; for i = 1,2ls Si^n = 1- The minimum value of;^ is then given by ^ = 0, giving 



(A.3) 
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Applying the chain rule results in 



( 4 



i=2 n \j=\ 



^Vl,n + ^bij(Vj^n - n,n) 
7=2 



7=1 



7=2 



where the last two sums over j can be evaluated as 

4 {A 



i=2 n Vj=1 



7=2 



{SikVln - Si^nSlk(vi,n " Vl,„)) = 0. 



(A.4) 



This expression splits naturally into two separate cases: 
Fory^= 1, / = 2,3,4: 
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andfory^> 1, /= 1,2,3,4 
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The solution of this system gives the least-squares fit to the data, which can be renormalized by choosing an appropriate value for A. The natural choice fo r A ii 



2, as the maximum intensity measured by one of the beams of the polarimeter is ll for unpolarized light. The one remaining element bn is constrained by Eq. A.3 



